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ABSTRACT 


The  simultaneous  effect  of  ion  heating  and  spatial  diffusion  due 
to  ICRF  heating  in  tokamak  geometry  is  investigated  from  the  single 
particle  viewpoint.  Transitions  between  passing  and  trapped  tra¬ 
jectories  are  taken  into  account  in  the  mapping  equations.  Specific 
results  are  given  for  typical  experimental  devices. 
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The  heating  of  plasma  ions  by  fast  Alfven  waves  of  frequency  w  in  the 

neighborhood  of  the  ion  cyclotron  frequency  (ICRF)  is  presently 

considered  to  be  an  attractive  scheme  for  achieving  thermonuclear  tempera- 

1  2 

tures  in  Tokamak  confinement  devices.  Previous  theoretical  studies  * 
of  this  topic  have  dealt  with  important  questions  related  to  the  efficiency 
of  the  heating  mechanisms  (fundamental,  harmonic,  minority),  and  the 
associated  evolution  of  the  ion  distribution  function  in  velocity 
space.  However,  such  studies  have  not  considered  the  intrinsic  diffusion 
in  configuration  space  which  must  always  accompany  the  velocity  space 
changes  produced  by  the  waves.  In  an  actual  experiment  if  the  spatial 
diffusion  rate  is  relatively  fast  compared  to  the  heating  rate,  then  the 
effect  of  the  heating  scheme  can  be  quite  deleterious  to  plasma  confine¬ 
ment  either  due  to  the  direct  escape  of  single  particles  (i.e.,  a  density 
loss)  or  because  of  the  indirect  recycling  of  wall  impurities.  The 

degrading  of  plasma  confinement  during  ICRF  heating  has  been  documented 

3  4  5 

experimentally  in  the  ST  ,  ATC  ,  and  Macrotor  Tokamaks.  Recent  experi¬ 
ments  in  PLT^  have  not  shown  any  deleterious  effects  of  the  RF  on  confine¬ 
ment,  as  expected,  because  of  the  large  current  and  size  of  this 
device  and  the  low  RF  power  levels  applied  thus  far. 

The  present  study  is  concerned  with  the  calculation  of  the  simultan¬ 
eous  change  in  the  energy  and  radial  position  produced  by  resonant 
heating  (w  -  k,,v„  =  nf2^,  n  =  1,2...)  of  single  ions  confined  in  a  Tokamak. 
Specifically,  this  work  yields  the  self-consistent  heating  and  spatial 
diffusion  rates  in  the  absence  of  collisions.  The  role  of  collisions 
will  be  treated  in  a  separate  st  'c'y.  An  important  feature  retained  in 
the  present  calculation,  and  generic  of  Tokamaks,  is  the  continuous 
transition  between  passing  and  trapped  particle  orbits  produced  by  the 
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preferential  increase  in  perpendicular  energy  due  to  cyclotron  absorption. 

The  physical  model  behind  the  present  calculation  hinges  on  the 
existence  of  a  finite  resonance  layer  inside  the  plasma  where  the  con¬ 
dition  b)  -  k„v„  =  n  ST  is  satisfied  for  a  wave  of  parallel  wavenumber  k„ 
and  a  parallel  particle  velocity  vlt.  This  arises  from  the  R“*  dependence 
of  the  toroidal  field,  i.e.,  n^R-1  where  R  is  the  major  radius.  For 

simplicity,  we  take  the  layer  to  pass  through  the  center  of  the  plasma 
minor  cross  section.  Because  c£  the  rotational  transform,  a  given  ion 
spends  only  a  finite  time  inside  the  layer.  During  this  time  the  ion 
magnetic  moment  y  experiences  a  net  change  Ay  due  to  the  RF.  After 
leaving  the  layer,  the  ion  is  no  longer  in  resonance  with  the  wave  and 
proceeds  to  follow  a  Tokamak  orbit  dependent  upon  the  single  particle 
invariants,  i.e.,  the  energy,  y,  and  the  toroidal  canonical  angular 
momentum  P  After  the  particle  circulates  half  way  around  the  magnetic 
surface,  it  encounters  the  resonance  layer  again.  At  this  point  the 
strength  of  the  RF  kick  received  by  the  particle  depends  on  the  phase 
difference  experienced  by  the  wave  while  the  particle  has  been  circulating. 
The  repeated  random  scattering  in  y  and  the  associated  orbit  changes  are 
responsible  for  the  simultaneous  heating  and  spatial  diffusion  of  the 
particle. 

Using  a  local  orthogonal  coordinate  system  (x„,  Xj_,  x3)  aligned 
with  the  total  magnetic  field  B,  the  Lorentz  force  equation  can  be 
averaged  over  the  particle  gyro-angle  to  yield  the  evolution  equations. 

y1/2  =  eE+  (2mB)_1/2  cos  (k^xj  cos  6  (1) 

fi  =  o)-  k„v„  -  fii  -  eE+  (2mBy)"1//2  cos  (k^x^)  sin  6  (2) 
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d 

dt 


v 


It 


-  - ^  §i 


(3) 


where  e,  m  refer  to  the  ion  charge,  and  mass,  and  E+  is  the  amplitude  of 
the  left-hand  polarized  component  (resonant)  of  the  wave  having  perpen¬ 
dicular  wave  number  k x;  6=6-  (k,,x„  -  ait)  where  6  is  the  gyroangle 
between  vx  and  the  normal  to  B. 

The  important  changes  described  by  Eqs.  (1)  -  (3)  occur  at  the 
crossing  of  the  resonant  layer,  where  abiabaticity  is  broken;  outside 
the  layer  these  equations  predict  small  oscillations  which  are  neglected 
in  this  calculation.  To  extract  the  jump  Ay  at  the  layer,  Eqs.  (1)  and 
(2)  are  combined  into  the  exact  expression 


d 

dt  ' 

y1^2exp[i (5- 

' 

dt'  (ai  -  k„v„  -  0))] 

=  eE. 


-1/2 

(2mB)  cos(kj_xx)exp[-i 


dt'(o)  - 


k„v„ 


-.0)] 


(4) 


Observing  that  the  left  side  of  Eq.  (4)  is  of  the  form  g(t)exp[if (t)] , 
where  g(t)  varies  slowly  in  time,  the  method  of  stationary  phase  can  be 
used  to  obtain  the  value  of  y  =  y(t)  after  passage  through  the  n-th 
harmonic  resonance. 


y(t)  =  |  y1/,2(o)  ±  eE+  (2mB) (At)exp 


i[±  n/4  -  <l>] 


(5) 


where.  At  =  cos(kxxQ)T  for  n  odd,  and  At  =  sin(kxxQ)T  for  n  even,  and 


T  =  Jn.1(kxp)(Tr/2)1/2  [(dfl/dxJCnv,,  -  ck„y/e)]-1/2 


*  = 


dt'  (to  -  k„v„  -  nfl) 


J  .00 


(6) 

(7) 
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J  is  the  Bessel  function  of  order  n,  x  the  guiding-center  position, 
n  o 

p  the  gyro  radius,  c  the  speed  of  light.  All  the  spatially  dependent 

terms,  B,  x  in  Eq.  (5),  p,v,„  u  in  Eq.  (6)  and  t  in  Eq.  (7)  are 

o 

evaluated  at  the  resonance  point.  In  Eq.  (5)  the  (+)  refers  to  n  odd, 

(-)  to  n  even.  Physically,  At  represents  the  effective  time  spent  by 
the  ion  inside  the  resonant  layer. 

To  calculate  the  phase  change  At/ *  between  resonances,  one  uses  the 
small  inverse  aspect  ratio  e  =  r/R  expansion  for  the  toroidal  field,  B^ 

=  B  (1  -  ecosa) . 

At/*  =  2(e(oRq/v„a)  (1+a)1^2  [E(a2/2,<)  -  E(a1/2,<)  -  (F(q2/2,k) 

-  F(a1/2,K))/(l+a)]  -  k()R  q  (a2  -  c^)  (8) 

where  a  =  2€Boy/mv„2,  k2  *  2a/(l+a),  q  is  the  safety  factor,  and  F(a,te), 

E(<x,*0  are  the  incomplete  elliptic  integral  of  the  1st  and  2nd  kind. 

In  Eq.  (8)  v„  is  evaluated  at  the  resonant  point  (a=a1).  A 
passing  particle  has  a <1  and  thus  <  <1.  The  appropriate  expression  of 
A tji  for  the  trapped  particle  can  be  found ^  from  the  continuation  formula 

of  the  elliptic  integral  from<>l  to  <  <1. 

Having  found  analytic  expression  for  the  changes  in  (y,'|<),  an 
appropriate  set  of  mapping  equations  can  be  found.  Because  in  Tokamak 
geometry  dB/dx„  changes  sign  upon  successive  passages  through  the  resonant 
layer,  two  different  sets  of  difference  equations  must  be  used.  In 
addition,  because  the  RF  can  induce  changes  between  passing  and  trapped 
orbits,  different  expressions  must  be  used  as  an  ion  crosses  the  trapping 


boundary. 
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The  change  of  the  ion  orbit  can  be  found  from  the  conservation 
of  the  total  energy,  p  and  P  between  resonances,  i.e.. 


2£pB  .. 

P,  =  omRv,,  (1  +  - j  cosa)  '  cosB  -  e  4>(r)/c 

mv„ 


(9) 


where  e*  ±  1  refer  to  co-flowing  and  counter-flowing  particles,  tan  g  =  B^/B^. 
From  Eqn.  (9) »  r  (a,p)  can  be  found  for  given  values  of  P^ ,  v„  and  the 
functional  form  of  the  poloidal  flux  function  Yfr).  After  updating  p-* 

P+Ap  using  Eqs.  (5),  (8) ,  the  orbit  change  can  be  obtained  from  Eq. 

(9). 


In  a  sample  numerical  calculation  one  considers  several  protons 
having  equal  p  initially  and  distributed  uniformly  over  the  phase  (0, 

2tt)  .  Then  Eqs.  (5),  (8),  (9)  are  iterated  by  a  computer  for  each 
particle  and  the  results  are  averaged  over  the  initial  phases.  For 
simplicity  a  cylinder  having  a  parabolic  current  profile  is  used  to 
obtain  t(r).  Fig.  (1)  shows  an  example  of  successive  changes  in  the 
averaged  total  energy,  while  Fig.  (2)  demonstrates  the  corresponding 
change  of  guiding  center  minor  radius  Ar  as  the  ion  passes  through  the 
equatorial  plane  (a  =  0) .  The  machine  parameters  correspond  to  the  UCLA 
Macrotor  for  3  different  RF  power  levels.  (The  relationship  between  E+ 
and  RF  power  can  be  found  with  the  aid  of  Ref.  (9).)  The  linear  dependence 
of  (Ar)  on  the  transit  number  seen  in  Fig.  (2)  demonstrates  that  the 
radial  motion  is  diffusive.  The  heating  rate  and  diffusion  coefficient 
can  be  obtained  from  the  slope  of  these  figures.  The  radially  averaged 
values  for  Macrotor  and  PLT  are  shown  in  Figs.  (3)  and  (4).  The  idealized 
parameters  used  to  generate  these  results  are  shown  in  Table  I.  All  the 
other  parameters  are  fixed  in  such  a  way  that  Fig.  (3)  and  (4)  represent 
first  harmonic  (n  =  1)  heating  of  protons  (minority)  in  a  deuterium  (majority) 
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plasma.  The  single  particle  confinement  time  and  maximum  attainable 

energy  can  be  found  by  examining  such  figures.  For  instance,  in  high 

current  PLT  at  300  kw  RF  power  it  takes  a  proton  45msec  to  attain  a 

perpendicular  energy  of  37.5  KeV  while  it  walks  out  radially  5cm. 

Similar  studies  have  been  made  for  second  harmonic  heating.  Since  this 

2 

process  is  proportional  to  (kxp)  ,  as  shown  in  Eq.  (6),  the  correspon¬ 
ding  heating  rate  and  diffusion  coefficient  are  smaller  than  the  funda¬ 
mental  case  at  least  by  an  order  of  magnitude. 

In  summary,  a  scheme  has  been  developed  which  permits  the  calcu¬ 
lation  of  the  simultaneous  ion  cyclotron  resonance  heating  rate  and  the 
associated  real  space  diffusion  coefficient  in  Tokamak  geometry.  The 
diffusive  nature  of  the  heating  and  radial  excursion  have  been  demon- 
strated,  and  its  magnitude  found  to  scale  as  |E+|  .  The  preferential 
transfer  of  perpendicular  energy  due  to  cyclotron  resonance  produces  an 
enhanced  banana  diffusion  rate  which  limits  the  maximum  energy  attainable 
by  single  particles.  Such  an  effect  can  be  quite  severe  for  small 
and/or  low  current  machines  and  should  be  included  in  sophisticated 
codes  that  attempt  to  explain  experimental  heating  result. 

The  authors  would  like  to  express  their  thanks  to  Professor  Y.  C.  Lee 
for  valuable  advice  on  this  problem. 
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Table  1.  Idealized  parameters  used  in  generating  Figs.  (3)  and  (4).  T„ 
is  the  initial  parallel  energy  at  the  resonance  layer. 


Machine 

Major 

Radius 

(cm) 

Minor 
Radius 
(cm)  _ 

B 

0 

(kG) 

I 

(KA)  _ 

Cavity 

Q 

T„ 

(eV) 

Macrotor 

100 

30 

2 

60 

30 

100 

PLT 

130 

40 

30 

600 

30 

500 

Fig.  1 


Fig.  2 


Fig.  3 


Fig.  4 


Figure  Captions 

Secular  increase  in  the  average  ion  energy  due  to  transits 
through  the  fundamental  resonance  (n  =  1)  for  different  RF 
power  levels.  Parameters  correspond  to  protons  in  the  Macrotor 
tokamak. 

Secular  change  in  the  guiding  center  minor  radius  Ar  as  the 
ion  passes  through  the  equatorial  plane.  This  radial  diffusion 
is  obtained  self- consistently  with  the  heating  shown  in  Fig.  1. 

Radial  and  phase  averaged  heating  rate  (solid  curve)  and 
spatial  diffusion  coefficient  (dashed  curve)  for  minority 
protons  in  Macrotor. 

Radial  and  phase  averaged  heating  rate  (solid  curve)  and 
spatial  diffusion  coefficient  (dashed  curve)  for  minority 


protons  in  PLT. 


Total  energy  (KeV) 


Fig.  1 


UJ3Z(JV) 


Rote  (eV/ms) 


Diffusion  Coefficient  (cm  /ms) 


Heating  Rate  (eV/ms) 


r- 

•4 

•» 

•4 

OO  '■«** 

O 

M 

X 

rH 

ai 

X 

DO 

c 

• 

ai 

(/) 

X 

V 

*-  ■ 

0-8 

M 

c 

o 

• 

•> 

r — n 

pH  • 

3 

V 

u. 

id 

o 

c/> 

Ui 

c 

o 

' — '  X 

id 

IL 

* 

4-» 

o 

• 

3 

•  .  T) 

•H 

oo 

w 

V 

Q 

•4 

>4 

M  D 

V) 

<- — « 

id 

3 

t — '  -H  •* 

X 

c 

ai 

O 

•  LL 

•4 

D 

M 

V  *H 

a* 

• 

Q 

C 

O  o  3  </> 

o 

r— < 

•d 

X 

• 

M 

V 

•d 

• 

XI  B 

• 

■4 

X 

*5 

00  00  ^  DO 

• 

</) 

x 

V  “ 

2 

O 

e 

3 

t — s 

6  Xi 

•J 

<7l 

V 

• 

<71  C7l  U.  DO 

H 

3 

3 

(L  c 

rH 

V 

c 

d 

v  O 

• 

rH 

X 

♦ 

2 

e 

*— *  r~~*  id 

•0 

X 

U. 

o 

• 

X 

LO 

X 

X 

>  V) 

CL 

' — / 

q 

< 

V 

W  t^,  X 

• 

O 

M 

/-V  V) 

id 

V 

ai 

O 

o 

• 

X 

o 

< 

id 

• 

<0  3 

H 

• 

a. 

•H 

Z  A 

3 

M 

B 

V 

XX  * 

• 

3 

2 

x  id 

•* 

< 

• 

' — v 

id 

C 

V 

•H 

V 

M 

M  M  <0  *J 

I 

2 

H 

u  o 

c 

• 

•  c 

-  B 

«d 

X 

Q 

•4 

id  id  u  l* 

« 

X 

•» 

id 

•  rH 

•4 

C  Id 

M 

C  *H 

*  £ 

c 

rH 

• 

3  3  -h  V 

3 

• 

V 

M 

M  • 

X 

c 

•  rH  CL 

V 

o  s 

A 

/ — '  V 

• 

o 

V 

cs  c  to  e 

rH 

3 

U. 

V 

*->  T. 

• 

c 

x  id 

43 

10 

X 

Ol  u 

O 

B 

c 

id  id  X  M 

V 

X 

V) 

• 

OS 

•d 

3 

O 

3  • 

D 

r~-  v 

*H 

e 

• 

*^3X3 

•  H 

3 

» 

o 

X  • 

f— 

M 

D 

id  ^ 

M 

cn  a 

a 

Q 

V 

CO 

a.  o 

LL 

C 

c 

id 

<  3 

u7 

H  » 

O 

Q 

3 

iH 

DO 

»  •» 

*5 

o 

E 

• 

3 

id 

o 

3 

*4 

-*  * 

C 

• 

X 

X  C  O  - 

O 

(0 

3 

DO  3 

< 

C 

X 

•  >» 

tn 

*  C 

V 

o 

C3 

• 

X 

>  id  x>  • 

•rH 

-> 

3 

X 

c  c 

V 

<  o 

* 

V 

s:  id 

• 

M  V 

32 

(L 

V  B  3 

M 

,3 

<d 

(J 

•h  id 

- 

X 

• 

DO 

t/> 

rH 

V) 

V  X 

*» 

X  r-^  -rt  C  • 
a>  v  ’H  ^ 

4-» 

4h 

o 

to 

M 

M 

V 

X 

*  id 

V 

id 

X 

• 

C 

• 

«/) 

O 

O 

V  * 

V 

3 

c  /, 

f-H 

M 

a 

fe 

e  - 

X 

o 

u 

id 

-J<  ^  3:  O 

V 

<0 

■ 

• 

V 

3 

r  i 

o 

cd 

o 

u!7  .,H 

•  S  u 

H 

& 

a>  l>  oo 

rH 

id 

2 

7* 

C  3 

V 

v>  - 

c  •  3 

O  O  -H 
V)  -r— * 

*  :  « 
«d  -  f— 

q  a 


a  u  .h  (i 

3  ..1.2 

-H  3  O 

a  a>  i/i  m 
M  M 

<D  C  -  0> 

4)>H£ 
X  4-»  id  D 

5t/>J  c 


'«  • 


do  in 
q 


M  •  X 
D  3  .q 
o 

<U  -  I 
*-h  £  D 
UJ  ffl  <u 
w  e  v> 


OS 

B 

U  r-H 

3 

A 

A  A  4> 

u  oc 

•  rH 

X 

X 

DO 

O 

»rH 

id  id 

X 

X 

rH 

3 

id 

(— 

0) 

>, 

r— H 

♦rH 

<U  rH 

O 

c 

C  H  .4 

V 

D 

jr 

V 

c 

3 

DO 

3c  V 

X 

• 

CL 

z 

rH 

X 

X 

X 

X 

V 

•n 

rH  X 

V 

o 

o  3  m 

•H  ?.  J-l 

x  • 

V 

3 

X 

•  rH 

O 

A 

t/> 

• 

♦  rH 

t/1 

o 

M 

• 

X 

cd 

c 

id 

x  id 

c 

•rH 

u 

q 

x 

M 

A 

rH 

t/1 

M  V 

t/l 

2 

D 

u 

a 

q 

O 

• 

f— 

3 

c 

H 

D 

V 

4-» 

4J  c  Q 

M 

DO 

S— H 

• 

O 

o 

o 

V 

•rH  0C 

M 

• 

rH 

•rH 

o 

D 

< 

id 

td 

V 

T3  vi 

DO 

id 

o  V 

U  X 

id 

D 

l/) 

Pd 

a) 

3 

V 

• 

3 

P 

D 

• 

M 

o 

A 

E 

p 

X 

• 

rH  C 

O 

> 

•ox- 

id  xi 

2 

D 

•rH 

D 

o 

rH 

£  r-H 

V) 

3 

S 

V) 

2T. 

D 

M 

• 

X 

V) 

3 

H 

o  ^ 

B 

5h 

f4  «  S 

C 

• 

c 

id 

H 

2 

DO  id 

id 

A 

O 

O 

O 

3 

id 

cC 

• 

u 

O 

V 

V  .  0) 

c 

V 

V 

3 

o 

X 

C  cj 

X 

A 

«d 

rH 

q 

A 

V 

rt 

O 

rH 

u 

>4 

c 

X 

to 

D  OS  X 

X  A 

a 

£ 

2 

3 

•d  *rH 

q 

6 

«5 

DO 

id 

rH 

A 

V 

X 

V 

4h 

id  O 

c 

X 

c 

X  M 

id 

•  rH 

A 

A 

rH 

U 

X  l/l 

q 

<d 

•«H 

3 

id 

£ 

X 

X 

M 

c 

3 

M 

t-H 

o 

HH  V)  A 

O 

X 

M 

z 

> 

DO 

V 

•rH 

X 

o 

MS 

•r- » 

■rH 

•rH 

•rH 

V 

M 

<D 

e 

<2 

•4 

V 

t*H  D 

•H  4-> 

2  M 

M 

Q) 

<D 

P 

c 

•H 

D 

C  X 

M 

id 

o 

Q 

**“t 

a> 

X 

id 

X 

o 

V 

O  O 

c 

V 

id  M  D 

id  m 

3 

S' 

M 

> 

o 

X 

V 

V  P,  D 

• 

X 

M 

id 

V 

D 

X 

to 

r-H 

•rH 

> 

B  M  V 

B  *H 

</) 

V 

c 

>  9 

u 

as 

O 

id 

P 

M 

• 

<u 

V 

M 

V 

V)  <J 

id 

V)  3  X 

V)  2 

X 

X 

M 

A 

DO 

•H  V 

V 

H- 

e 

X 

c 

• 

X. 

fj 

X 

o 

X 

<7  X 

X 

32 

id  u  u 

id 

X 

Pd 

O 

A 

►— < 

id 

M  O 

f-H 

• 

H 

t^ 

• 

•  rH 

X 

c 

D 

> 

•H  CJ 

D 

rH  D 

r-H  U 

X 

M 

X 

o 

2 

Q 

X 

3 

id 

id 

f— 

X 

• 

a> 

•H 

• 

B 

•H 

id 

Ou  X-H 

CL  -rH 

O 

o 

•d 

D 

• 

l  X 

A 

D 

rH 

D 

V) 

X 

3 

$H 

C 

c 

2 

•d  C 

rH 

B 

l  X  u 

D 

3 

M 

D 

V) 

X 

M 

u 

e  o 

V 

A 

rH 

V 

X 

DO 

q 

X 

<D 

•H 

C3 

• 

£5 

•rH 

to 

h  CL 

B  V 

V 

M 

V 

M 

3 

id 

•  rH 

id 

a> 

V 

V 

CO 

V) 

q 

o 

X 

X 

3 

X 

id 

v  r 

O  C 

M 

•rH 

1) 

O 

CO 

CL 

M 

a>  rH 

M 

q 

OJ 

*-> 

V 

•  rH 

r 

X 

c 

B 

3 

Q  * 

rH 

V)  A  . 

M  DO 

•  H 

2 

D 

D 

x  id 

X 

3 

X 

M 

<D 

X 

D 

id 

o 

id 

•rH 

c 

C 

4-» 

CL 

id  C  0u 

X  jA 

rH 

id 

U 

A 

1 

U 

c 

•d 

M 

DO 

3 

id 

B 

r 

CJ 

<D 

3 

•4 

DO  V 

to 

x  o 

£ 

•rH 

u 

2 

V 

M 

o; 

X  M 

6 

2£ 

CO 

•rH 

D 

U 

V 

3 

t/i 

in 

CO 

r-H 

on 

e  > 

c 

B 

•H  T 

D 

X 

•rH 

D 

V 

V) 

f-H 

O  3 

o 

X 

id 

o 

*-» 

<d 

D 

rH 

u. 

C 

•H  <fH 

►H 

o 

X  D  a 

X  id 

id 

M 

f-H 

V 

X 

D 

X 

O 

M 

• 

• 

M 

X 

•rH 

O 

r-H 

C 

id 

c 

o 

Ph  rH 

M 

o  <d  v 

DO 

4-» 

D 

Q 

o 

c 

0)  3 

X 

CO 

X 

A 

V 

c 

3 

X 

•  rH 

c 

o 

to 

u 

P-< 

c 

<4H 

3  3 

•rH  C 

CO 

V 

c 

V 

>4 

m 

X 

CL 

V 

•  iH 

rH 

o 

o 

M 

3 

id 

o 

0  7  0 

X  .rH 

E 

o 

V 

B 

B 

M 

CL  O 

q 

• 

D 

6 

> 

X 

X 

X 

X 

•  rH 

D 

O 

o 

M  V 

M 

rH 

e  X  jr 

V) 

P 

•rH 

3 

3 

•  H 

Id  D 

o 

3 

A 

•  rH 

<1) 

•d 

M 

1 

•M 

O 

V 

a> 

h-  X 

D 

•H 

o 

3  X 

M 

V) 

•rH 

X 

M 

D 

rH 

•rH 

q 

to 

f— 

id 

o 

o 

X 

<d 

<U 

c 

r— H 

D 

O 

id 

•H  C  DO 

1  D 

O 

• 

<n 

c 

D 

o 

V 

•rH 

rH  T* 

D 

* 

B 

q 

•rH 

D 

q 

V 

X 

V 

rH 

D 

DOC 

O  *rH 

M 

/ - 4 

•  rH 

V 

o 

CO 

Pd  rH 

o  v 

id 

id 

3 

V 

X 

M 

X 

D 

id 

DO 

|h 

in 

DO 

X  3 

o 

o 

Id  »H  *rH 

rH  U 

M 

Cl 

K 

B 

D 

o 

O  D 

•H 

B 

rH 

o 

DO 

•rH 

o 

id 

M 

c 

a> 

o 

X  D  C 

X 

D 

rH  V)  M 

•  rH 

•rH 

X 

< 

•rH 

o 

X  X 

co 

D 

3 

•  rH 

o 

• 

• 

•  rH 

3 

D 

O 

•H 

*~> 

C 

e 

id 

id 

CJ 

V 

3  3  id 

X  D 

2 

Oi 

a 

X 

X 

rH  *rH 

•d 

•' — » 

LJ 

D 

X 

4—-V 

X 

B 

>4 

CO 

X 

X 

id 

o 

3 

C 

• 

c 

B  X  0) 

O  V) 

rH 

V 

i 

CL 

c 

V 

•2J 

as 

id 

D 

DOO 

DO 

4-1 

O 

id 

M 

» — < 

X 

id 

O  DO 

c 

/ - 4 

DO 

.H  X  H 

id 

o 

DO 

V 

z 

o 

V) 

H 

id 

V 

M 

00 

q 

q 

•rH 

M 

X 

x  < 

p 

c 

o 

Cl 

t/)  »H 

M  X 

• 

•rH 

M 

V 

B 

Cl 

•  H 

c 

D  3 

t> 

e  -j 

V 

Ci 

o 

« 

rH 

D 

> 

a> 

X 

3 

3  -h 

►H 

Q  V 

a>  o 

V) 

M 

M 

id 

M 

3 

X 

D 

V 

id  co 

•  rH 

• 

</) 

q 

rH 

X 

*rH 

U 

id 

l 

x: 

V 

OS 

rH  *J 

1 

Cl 

rH  rH 

D  O 

3 

D 

V 

X 

iS 

3 

• 

<71 

O 

D 

Pd 

4-> 

H 

rt 

• 

X 

v — ' 

V 

X 

<D 

V 

+j 

3 

X 

v  id 

c 

rH 

c 

Id  D  x 

D  D 

•rH 

V 

X 

u 

f-H 

V 

c 

t/) 

V 

rH 

> 

o 

X 

id 

rH 

B 

E 

V 

D 

C 

•rH  V 

V 

4 - ' 

id 

D  C  3 

id  to 

3 

g 

• 

KJ 

O 

c 

HH 

• 

q 

2 

X 

V 

X 

X 

q 

o 

D 

X 

o 

•  rH 

X 

o 

o 

tL  X 

> 

•H 

C  V  o 

u 

rH 

E 

V 

> 

00 

V) 

c 

3  X 

DO 

A 

DC 

DO 

M 

V 

CO 

M 

f~ 

o 

c 

<D 

Mh 

M 

> 

V  rH  o 

IL 

X 

>  u 

<7» 

X 

cn 

D 

o 

X 

C  M 

id 

M 

*3 

•H 

id 

M 

q 

q 

u7 

X 

V 

O 

c 

u  <u 

r-H 

V 

o 

B  3 

V) 

o 

=3 

X 

4-> 

3 

rH 

rH 

u 

X 

id  id 

B 

V 

D 

• 

X 

3 

V 

o 

*-H 

X 

t/1 

•rH 

o 

•H  U 

-t 

X 

3 

•H  X  V 

o  u  x 

■H 

Cl 

V 

w 

3 

V 

u 

O 

X 

V 

to 

M 

X 

*H 

U 

DO 

o 

DO 

M 

V) 

•H 

u  c 

B 

MMX 

•H  a 

o 

X 

V 

rH 

X 

M 

to 

oc 

c 

DO  to 

M 

id 

CO 

xn-< 

X 

q 

D 

V 

•H 

C 

c 

V 

c 

D 

d  id 

V 

55 

V 

V  3  D 

U  13 

< 

X 

v — ' 

o 

V 

X 

V 

o 

c 

D 

X 

1 

X 

o 

V 

o 

3 

D 

D 

•  rH 

to 

•  rH 

P< 

o 

id 

u  c 

a*  A 

X 

> 

A 

X 

u 

L. 

44  C 

V) 

A 

KA 

D 

t  A 

D 

M 

OS 

V 

r* 

•H 

u 

•H  3 
i,  r* 

O 

At 

X 

KA 

Pd 

X 

u 

4-> 

3 

V 

r- 

M 

c 

3 

U1  V13  H 

*  jg  uj 

DO  43  OC 
05  M  U 
D  3  MQ 

a  v  m 
5-i  M  rt 
a>  a)  a 

>  CL  X 
•H  *3  0 
V  X  0) 

O  D  U  D 

5  «H  C  O 
H  If)  14  V 
^  p  M  X 
O  «  C  'W 
u  O  UJ  UJ 


cn  z: 

v  X  X 
ol  o  *  o  • 

Ifl  r-N 
C  CD  4)0) 


O  id  x  u  ' — ' 

r-H  M  3 
V  »H  •  M  L 
Xh  (fl  v  D 
U  -H  X  10  JO 
X  »C  B 
c  «a  «)  v 
O  D  ,c  > 
w  w  O  h  O 
2  2  D  S  Z 


V)  14  V 
•H  »H  *3 
t/>  O  -H 

X  * 

h  u5  Tj 

«5  C 

e  c  * 

<  *H  ffl 

MOV 
«S  D 

V  «-H  -H 


V)  4->  3 

x  id  o 

Xtf>- 
X  I 

m  • 

O  O  CL 

*  t  . 
Iso 


3  V  C 
♦H  JO  o 

b  -h 

<0  V  D 
d  u  id 

H  V  ^ 

3  a  3 

w  e 

V  -H 
OS  •  CO 

DO 

V  c 

e  3  « 
VGA 


«  -o  X  O  O  D 
X  41  rt  Q,-r4 
M  t<  H  K  W  U 


•HOC  «^H4JO  OO 

I  u  (i  in  in  m  in  ro 

e  a  o  e  .o  •  t: 

hjtj  4>  ♦->  —t  o  <  /-. 

o  c  «)  o.  o  o  <h 

E4J-HXOOC00O 
o  W  14-1  f—  DC  RJ  Ch 
M  -H  C  •  —  C 

IH  e  O  <  V)  UJ  41  V. •  o 

43  u  a.  O  VI 
,  VI  3  *  <  (X  X  -H 

W  W  41  M  3  X  Ih  (4 

o  -HC-a4CM«u*j 

HG£0l<ll|03a 


4NA9IOH 
(O  V)  m  IG  ^  t 
4t  tf  <J  >t  f  ^ 
I  I  I  I  I  I 

ouuouu 
Si  SI  c*  c-  &•  <x 

a.,  a.  a.  a.  a.  a. 


to 

q  -*-h 

c 

a 

u 

X 

u 

X 

o 

c 

rH 

3 

3 

DO 

M 

M 

rH 

£ 

3 

X 

c 

3  .O 

cd 

X 

•rH 

V 

V 

c 

R 

00 

£ 

•H  H 

•rH 

X 

id 

V 

rH 

« 

•H 

O 

Id 

•d 

id 

r-H 

c 

C 

1 

rO 

q 

D 

M 

fS 

-p 

o 

<71 

< 

X  CO 

X 

3 

X 

• 

es 

• 

u 

3 

X 

• 

Q 

X 

X 

X 

id 

> 

V 

X 

o 

tn 

• 

00 

V 

t.J 

rH 

t.J 

-  * 

= 

to 

1 

-1 

X 

S 

' — f 

1 

= 

= 

o 

: 

-3 

X 

" 

H 

r 

*r. 

= 

X 

W 

in  x 

X 

00 

Cl 

o 

rH 

ra 

t 

m 

X 

r- 

00 

Cl 

o 

r-H 

IM 

m 

-d- 

rf 

m 

in 

in 

in 

in 

m 

in 

m 

m 

m 

X 

X 

X 

X 

X 

rr 

Tf 

•d* 

t 

rf 

*± 

-r 

-t 

rt 

•d- 

-* 

a 

o  u 

2 

o 

a 

a 

O 

O 

<J3 

o 

o 

o 

a 

u 

a 

o 

a 

a 

a 

o 

CL 

X  X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X  X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Cl 

Cl 


£  m  ->  </l 

«  6  A 

"1  X>  rt  (%| 

CO  M  C 

r*  -  o  ■  jc. 

G-  Bl-rt  M  4J 

2  .  6  U1  h  -H 

g  **»  “1-H  II)  J 

"J  {?  <a  — ' 

WO  i — I  H  fj  (/) 

r1  *-*  g.  o  o>  e 

G.  CO  O  *-J  cO 

4->  CO  o 

*  O  n]  j(  cQ 

»-H  X)  c  S 

*-*  W  *H  C  .r-«  G 

-  J  SS 

G  w  £  w  a; 

w  oo-o  o 

Q-  C  *«H  I— «  •  *H 

’rt  *-J  O  (/)  *J 

n  M  •  >h  U?  j~'  oo 

£  JO  '->  »  £  >H 

G  G  o  cnu,  (u 

5  Q  «  cOQ£  O  c 

e  cr>  -4  rt  tu 

■r<  c  <-«  »rt 

t!  o  ^  CNQ  -Q  <4-. 

®  ►*>  „  O  <u  o 

8*  „  S  ^ 

,n  s  b,  o  <->  <o 

M  2  «  V<rt  rt  e 

o  3  o«  e  co 

g  r;  h  co  /)  o 
0  0.0  00  Si  3  s 

r!  £j  c  <->  in 

4-*  U.  -H  4)  >. 

O  C  rt  e  ja 

VC.  •>  44  4)  . 

C  MB,  O  </l 

C  -O  <U  *J  oo  rt 

O  C  n  cd  </)  CTi  B 

y  •  C  O  3  rt  <n 

“  *)  «  O  CO 

C  P  f— ♦  at  r-H 

O  ecu  x  o, 

to  *H  B  O 

fZ  4-t  .301^^ 

to  uy  m  cd 

-a  rt  -.-<<+-1  .  ^ 

*~4  44  c  rt  O  44  O 

44  o  CO  CO  rt  t— 

•74  O  B  C  N 

G-  <  "-1  T)  O  C  >*4 

*  ®-<-c  V  O 

<j  4)  xc  *->  rt  ue 

rt  rt  04  CO  CO  O  00 
rt  U  u  rt  >  (0  c 

St!  2  p  t  rt 

P.  {-4  .  •  £  4)  4-> 

00  *4  X.  rt  to  .  CO 

•S  J?  .  *^:S  «  ® 
f  ?•  £*  S°P  *  ? 

>g  ^  00  O)  O 

\©  >0  -0\0  r- 

"»  "»  4t4f  « 

I  •  Cl  I 

£g  gg  g 

O-  G.  0.0.  o. 


